Works on an information geometrodynamical approach to chaos 



Carlo CafarcB 

Department of Physics, State University of New York at Albany-SUNY, 
UOO Washington Avenue, Albany, NY 12222, USA 

In this paper, I propose a theoretical information-geometric framework suitable to characterize 
chaotic dynamical behavior of arbitrary complex systems on curved statistical manifolds. Specif- 
ically, I present an information-geometric analogue of the Zurek-Paz quantum chaos criterion of 
linear entropy growth and an information-geometric characterization of regular and chaotic quan- 
tum energy level statistics. 
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^ i I. INTRODUCTION 

. The study of complexity [l[ has created a new set of ideas on how very simple systems may give rise to very complex 
CNJ ' behaviors. Moreover, in many cases, the "laws of complexity" have been found to hold universally, caring not at all 
for the details of the system's constituents. Chaotic behavior is a particular case of complex behavior and it will be 
f 1 the object of the present work. 

In this paper I make use of the so-called Entropic Dynamics (ED) @. ED is a theoretical framework that arises 
from the combination of inductive inference (Maximum Entropy Methods (ME), Q) and Information Geometry 
(IG) [4(]. The most intriguing question being pursued in ED stems from the possibility of deriving dynamics from 
d . purely entropic arguments. This is clearly valuable in circumstances where microscopic dynamics may be too far 
removed from the phenomena of interest, such as in complex biological or ecological systems, or where it may just 
1— 1, be unknown or perhaps even nonexistent, as in economics. It has already been shown that entropic arguments do 
account for a substantial part of the formalism of quantum mechanics, a theory that is presumably fundamental [f| • 
V/ Perhaps the fundamental theories of physics are not so fundamental; they may just be consistent, objective ways 
to manipulate information. Following this line of thought, I extend the applicability of ED to temporally-complex 
^f*) j (chaotic) dynamical systems on curved statistical manifolds and identify relevant measures of chaoticity of such an 
\^ . information geometrodynamical approach to chaos (IGAC) . 

The layout of the paper is as follows. In the next Section, I give an introduction to the main features of our IGAC. 
In Section III, I apply my theoretical construct to three complex systems. First, I study the chaotic behavior of 
— 1 ■ an ED Gaussian model describing an arbitrary system of I degrees of freedom and show that the hyperbolicity of 
the non-maximally symmetric 2Z-dimensional statistical manifold M s underlying such ED Gaussian model leads to 
linear information geometrodynamical entropy (IGE) growth and to exponential divergence of the Jacobi vector field 
intensity. An information-geometric analogue of the Zurek-Paz quantum chaos criterion of linear entropy growth and 
an information-geometric characterization of regular and chaotic quantum energy level statistics are presented. 

I emphasize that I have omitted technical details that will appear elsewhere. However, some applications of my 
IGAC to low dimensional chaotic systems can be found in my previous articles [g, [D, l8|, [9( . Finally, in Section IV I 
present my conclusions and suggest further research directions. 
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II. THE INFORMATION GEOMETRODYNAMICAL APPROACH TO CHAOS: GENERAL 

FORMALISM 

The IGAC is an application of ED to complex systems of arbitrary nature. ED is a form of information-constrained 
dynamics built on curved statistical manifolds Ms where elements of the manifold are probability distributions 
{P(A|0)} that are in a one-to-one relation with a suitable set of macroscopic statistical variables {9} that provide 
a convenient parametrization of points on Ms- The set {&} is called the parameter space T>q of the system. 

In what follows, I schematically outline the main points underlying the construction of an arbitrary form of entropic 
dynamics. First, the microstates of the system under investigation must be defined. For the sake of simplicity, I 
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assume the system is characterized by an /-dimensional microspace with microstates {x k } where k — 1,..., I. These 
microstates can be interacting, non-interacting, or I may have no relevant information concerning their microscopic 
dynamics. Indeed, the main goal of an ED model is that of inferring "macroscopic predictions" in the absence of 
detailed knowledge of the microscopic nature of arbitrary complex systems. Once the microstates have been defined, 
I have to select the relevant information about such microstates. In other words, I have to select the macrospace of 
the system. For the sake of the argument, I assume that our microstates are Gaussian-distributed. They are defined 
by 2/-information constraints, for example their expectation values fi k and variances a k . 

i 

(x k ) = n k and (Ux k - (x k )) 2 )) 2 = a k . (1) 



In addition to information constraints, each Gaussian distribution p k (x k \fi k , o~ k ) of each microstate Xk must satisfy 
the usual normalization conditions, 

+oo 

J dx k p k {xk\n k , o-fe) = 1 (2) 

— oo 

where 

Pk (x k \fx k , a k ) = (2™ 2 fe p exp (- tZhZjjj^ j . (3) 

Once the microstates have been defined and the relevant (linear or nonlinear) information constraints selected, I am 

i 

left with a set of probability distributions p(X\Q) = J^J p k (x k \fi k> cr k ) encoding the relevant available information 

fc=i 

about the system where X is the /-dimensional microscopic vector with components (x\,...,xi) and O is the 21- 
dimensional macroscopic vector with coordinates (/x lv .., Hi\ <Ji,..., a{). The set {0} define the 2Z-dimensional space 
of macrostates of the system, the statistical manifold Ms- A measure of distinguishability among macrostates is 
obtained by assigning a probability distribution P (X|0) 9 Ms to each macrostate O . Assignment of a probability 
distribution to each state endows Ms with a metric structure. Specifically, the Fisher-Rao information metric g^u (O) 

1. 

3^(6) = J dXp(X\Q)d^ogp(X\e)d u \ogp(X\e), (4) 

with fj,,i> = 1,..., 21 and — , defines a measure of distinguishability among macrostates on Ms- The statistical 
manifold Ms, 

Ms = \p(X\Q) = l[ Pk (x k \ l x k ,a k )\, (5) 



fc=i 



is defined as the set of probabilities {p (X\Q)} described above where X G R 3JV , O € Vq = [2" M x X a ] 3N . The parameter 
space T>q (homeomorphic to Ms) is the direct product of the parameter subspaces X M and where (unless specified 
otherwise) X M = (-co, + oo)^ and X a = (0, +00)^. Once Ms and T>q are defined, the ED formalism provides 
the tools to explore dynamics driven on Ms by entropic arguments. Specifically, given a known initial macrostate 
Q(initiai) (probability distribution), and that the system evolves to a final known macrostate 0( final ), the possible 
trajectories of the system are examined in the ED approach using ME methods. 

I emphasize ED can be derived from a standard principle of least action (of Jacobi type) . The geodesic equations 
for the macrovariables of the Gaussian ED model are given by nonlinear second order coupled ordinary differential 
equations, 

^ + r ^wT = °- (6) 

The geodesic equations in ([6]) describe a reversible dynamics whose solution is the trajectory between an initial 0( inltlal ) 
and a final macrostate ©( fina1 '. The trajectory can be equally well traversed in both directions. Given the Fisher-Rao 
information metric, I can apply standard methods of Riemannian differential geometry to study the information- 
geometric structure of the manifold Ms underlying the entropic dynamics. Connection coefficients T^, Ricci tensor 
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R^v, Riemannian curvature tensor R^ vpiJ , sectional curvatures ICm S i scalar curvature 7Zm S i Weyl anisotropy tensor 
W^prj, Killing fields £ M and Jacobi fields J M can be calculated in the usual way. 

To characterize the chaotic behavior of complex entropic dynamical systems, I are mainly concerned with the signs of 
the scalar and sectional curvatures of Ads, the asymptotic behavior of Jacobi fields J M on A4s, the existence of Killing 
vectors £^ (or existence of a non-vanishing Weyl anisotropy tensor, the anisotropy of the manifold underlying system 
dynamics plays a significant role in the mechanism of instability) and the asymptotic behavior of the information- 
geometrodynamical entropy (IGE) Sm s (see ©)■ It is crucial to observe that true chaos is identified by the occurrence 
of two features: 1) strong dependence on initial conditions and exponential divergence of the Jacobi vector field 
intensity, i.e., stretching of dynamical trajectories; 2) compactness of the configuration space manifold, i.e., folding of 
dynamical trajectories. The negativity of the Ricci scalar TZmsi 

Km s = R^pag f " p g l/a = ^2 ( e P' e<J "> ' ( 7 ) 

implies the existence of expanding directions in the configuration space manifold Ai s . Indeed, since TZ-Ms is the sum 
of all sectional curvatures of planes spanned by pairs of orthonormal basis elements {e p = <9e p }, the negativity of the 
Ricci scalar is only a sufficient (not necessary) condition for local instability of geodesic flow. For this reason, the 
negativity of the scalar provides a strong criterion of local instability. Scenarios may arise where negative sectional 
curvatures are present, but the positive ones could prevail in the sum so that the Ricci scalar is non-negative despite 
the instability in the flow in those directions. Consequently, the signs of ICm s are °f primary significance for the 
proper characterization of chaos. 

A powerful mathematical tool to investigate the stability or instability of a geodesic flow is the Jacobi-Levi-Civita 
equation (JLC equation) for geodesic spread, 

D 2 J» „ d<d u n d<d a 

The JLC-equation covariantly describes how nearby geodesies locally scatter and relates the stability or instability of a 
geodesic flow with curvature properties of the ambient manifold. Finally, the asymptotic regime of diffusive evolution 
describing the possible exponential increase of average volume elements on M. s provides another useful indicator of 
dynamical chaoticity. The exponential instability characteristic of chaos forces the system to rapidly explore large 
areas (volumes) of the statistical manifold. It is interesting to note that this asymptotic behavior appears also in the 
conventional description of quantum chaos where the entropy (von Neumann) increases linearly at a rate determined 
by the Lyapunov exponents. The linear increase of entropy as a quantum chaos criterion was introduced by Zurek 
and Paz [10(. In my information-geometric approach a relevant quantity that can be useful to study the degree of 
instability characterizing ED models is the information geometrodynamical entropy (IGE) defined as [7| , 

S Ms (r) d = ^lim o logV A1s with V Ms (r) = - J dr' I J yfgd 2l <d J (9) 

o \m s / 

and g — |det (<? M „)|. IGE is the asymptotic limit of the natural logarithm of the statistical weight defined on Ai s and 
represents a measure of temporal complexity of chaotic dynamical systems whose dynamics is underlined by a curved 
statistical manifold. In conventional approaches to chaos, the notion of entropy is introduced, in both classical and 



quantum physics, as the missing information about the systems fine-grained state For a classical system, suppose 
that the phase space is partitioned into very fine-grained cells of uniform volume Av, labelled by an index j. If one 
does not know which cell the system occupies, one assigns probabilities pj to the various cells; equivalently, in the 
limit of infinitesimal cells, one can use a phase-space density p(Xj) — 4^-. Then, in a classical chaotic evolution, the 
asymptotic expression of the information needed to characterize a particular coarse-grained trajectory out to time r 
is given by the Shannon information entropy (measured in bits), 

S£ZS = -JdX P (X)lo g2 (p(X)Av) 

= -^2pjlog 2 pj ~ K.t. (10) 

3 

where p (X) is the phase-space density and pj — is the probability for the corresponding coarse-grained trajectory. 

^ciassica? ^ s the missing information about which fine-grained cell the system occupies. The quantity IC represents the 
linear rate of information increase and it is called the Kolmogorov-Sinai entropy (or metric entropy) (IC is the sum of 
positive Lyapunov exponents, IC — Xj ). IC quantifies the degree of classical chaos. 
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III. THE INFORMATION GEOMETRODYNAMICAL APPROACH TO CHAOS: APPLICATIONS 

In this Section I present three applications of the IGAC. First, I study the chaotic behavior of an ED Gaussian 
model describing an arbitrary system of I degrees of freedom and show that the hyperbolicity of the non-maximally 
symmetric 2Z-dimensional statistical manifold A4 S underlying such ED Gaussian model leads to linear information 
geometrodynamical entropy (IGE) growth and to exponential divergence of the Jacobi vector field intensity. I also 
present an information-geometric analogue of the Zurek-Paz quantum chaos criterion of linear entropy growth. This 
analogy is presented by studying the information geometrodynamics of ensemble of random frequency macroscopic 
inverted harmonic oscillators. Finally, I apply the IGAC to study the entropic dynamics on curved statistical manifolds 
induced by classical probability distributions of common use in the study of regular and chaotic quantum energy level 
statistics. In doing so, I suggest an information-geometric characterization of regular and chaotic quantum energy 
level statistics. 

As I said in the Introduction, I have omitted technical details that will appear elsewhere. However, my previous 
works (especially ([7|)) may be very useful references in order to clarify the following applications. 

A. Chaotic behavior of an entropic dynamical Gaussian model 

As a first example, I apply my IGAC to study the dynamics of a system with I degrees of freedom, each one described 
by two pieces of relevant information, its mean expected value and its variance (Gaussian statistical macrostates) . 
The line element ds 2 — g^ LV (6) dQ^d® u on M s is defined by, 

ds 2 = (—d/j, 2 . -\ — j^fe ) ' v ~ I'"-' ^- (11) 

k=i ^ °* ' 

This leads to consider an ED model on a non-maximally symmetric 2Z-dimensional statistical manifold M. s . It is shown 
that M s possesses a constant negative Ricci curvature that is proportional to the number of degrees of freedom of 
the system, Rm s = —I- It is shown that the system explores statistical volume elements on M s at an exponential 
rate. The information geometrodynamical entropy Sm 3 increases linearly in time (statistical evolution parameter) 
and is moreover, proportional to the number of degrees of freedom of the system, Sm s r ^°° IXt. The parameter A 
characterizes the family of probability distributions on M. s . The asymptotic linear information-geometrodynamical 
entropy growth may be considered the information-geometric analogue of the von Neumann entropy growth introduced 
by Zurek-Paz, a quantum feature of chaos. The geodesies on M s are hyperbolic trajectories. Using the Jacobi-Levi- 
Civita (JLC) equation for geodesic spread, I show that the Jacobi vector field intensity Jm 3 diverges exponentially and 
is proportional to the number of degrees of freedom of the system, Jm s t ^ j °° Zexp (Ar). The exponential divergence 
of the Jacobi vector field intensity J m s is a classical feature of chaos. Therefore, we conclude that 

1Zm s = —I, Jm 3 T ~°° Zexp(Ar) , Sm s t ^°° IXt. (12) 

Thus, TZm 3 , $M a an d Jm 3 behave as proper indicators of chaoticity and are proportional to the number of Gaussian- 
distributed microstates of the system. This proportionality, even though proven in a very special case, leads to 
conclude there may be a substantial link among these information-geometric indicators of chaoticity. 

B. Ensemble of random frequency macroscopic inverted harmonic oscillators 

In our second example, I employ ED and "Newtonian Entropic Dynamics" (NED) 9]. In this NED, we explore 
the possibility of using well established principles of inference to derive Newtonian dynamics from relevant prior 
information codified into an appropriate statistical manifold. The basic assumption is that there is an irreducible 
uncertainty in the location of particles so that the state of a particle is defined by a probability distribution. The 
corresponding configuration space is a statistical manifold M s the geometry of which is defined by the Fisher- Rao 
information metric. The trajectory follows from a principle of inference, the method of Maximum Entropy. There is 
no need for additional "physical" postulates such as an action principle or equation of motion, nor for the concept 
of mass, momentum and of phase space, not even the notion of time. The resulting "entropic" dynamics reproduces 
Newton's mechanics for any number of particles interacting among themselves and with external fields. Both the 
mass of the particles and their interactions are explained as a consequence of the underlying statistical manifold. 
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In my special application, I consider a manifold with a line element ds 2 = g^ (9) dQ^dQ 1 ' (with fx, v = 1,..., V) 
given by, 



ds 2 = [i - $ (6)] v (9) de^e", $ (6) - Y. u x ( 13 ) 

where 



fe=i 



«fc(flfc) = -5^,e fc =fl fc ( a ). (14) 

The geodesic equations for the macrovariables (s) are strongly nonlinear and their integration is not trivial. How- 
ever, upon a suitable change of the afhne parameter s used in the geodesic equations, I may simplify the differential 
equations for the macroscopic variables parametrizing points on the manifold M. s with metric tensor g^. Recalling 
that the notion of chaos is observer-dependent and upon changing the affine parameter from s to r in such a way that 
ds 2 =2(1- $) dr 2 , I obtain new geodesic equations describing a set of macroscopic inverted harmonic oscillators 
(IHOs). In order to ensure the compactification of the parameter space of the system (and therefore A4 S itself), we 
choose a Gaussian distributed frequency spectrum for the IHOs. Thus, with this choice of frequency spectrum, the 
folding mechanism required for true chaos is restored in a statistical (averaging over to and r) sense. Upon integrating 
these differential equations, I obtain the expression for the asymptotic behavior of the IGE Sm s i namely 

i 

S M . (t) T ~°° At, A = (15) 
»=i 

This result may be considered the information-geometric analogue of the Zurek-Paz model used to investigate the 
implications of decoherence for quantum chaos. In their work, Zurek and Paz considered a chaotic system, a single 
unstable harmonic oscillator characterized by a potential V (x) = — ^ (fl is the Lyapunov exponent), coupled to 
an external environment. In the reversible classical limit [12J, the von Neumann entropy of such a system increases 
linearly at a rate determined by the Lyapunov exponent, 

fijSSSW^Or, (16) 
with Q playing the role of the Lyapunov exponent. 



C. Information geometrodynamics of regular and chaotic quantum spin chains 

In my final example, I use my IGAC to study the entropic dynamics on curved statistical manifolds induced by 
classical probability distributions of common use in the study of regular and chaotic quantum energy level statistics. 
Recall that the theory of quantum chaos (quantum mechanics of systems whose classical dynamics are chaotic) is not 
primarily related to few-body physics. Indeed, in real physical systems such as many-electron atoms and heavy nuclei, 
the origin of complex behavior is the very strong interaction among many particles. To deal with such systems, a 
famous statistical approach has been developed which is based upon the Random Matrix Theory (RMT). The main 
idea of this approach is to neglect the detailed description of the motion and to treat these systems statistically 
bearing in mind that the interaction among particles is so complex and strong that generic properties are expected to 
emerge. Once again, this is exactly the philosophy underlining the ED approach to complex dynamics. It is known 
[l3| that integrable and chaotic quantum antiferromagnetic Ising chains are characterized by asymptotic logarithmic 
and linear growths of their operator space entanglement entropies, respectively. In this last example, I consider the 
information-geometrodynamics of a Poisson distribution coupled to an Exponential bath (spin chain in a transverse 
magnetic field, regular case) and that of a Wigner-Dyson distribution coupled to a Gaussian bath (spin chain in a tilted 
magnetic field, chaotic case). Remarkably, I show that in the former case the IGE exhibits asymptotic logarithmic 
growth while in the latter case the IGE exhibits asymptotic linear growth. In the regular case, the line element 
^Ltcgrabic = d 4oisso„ + d 4x P onc„tiai on the statistical manifold M s is given by 

^ s intcgrablc = ^ 2~°^S (1?) 

I 1 A Ms 

where the macrovariable \x A is the average spacing of the energy levels and \i B is the average intensity of the magnetic 
energy arising from the interaction of the transverse magnetic field with the spin ^ particle magnetic moment. In 
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such a case, I show that the asymptotic behavior of i 5^ ltc s iablc ) j s sub-linear in r (logarithmic IGE growth), 

5 Mcgrablc) (r) r^oo logr _ 

Finally, in the chaotic case, the line element dSchaotic = ^ s wigncr-Dyson + <^ s Gaussian on the statistical manifold A4 S is 
given by 

4 12 
^chaotic = — d ^'l + —d^B + — d<J B ( 19 ) 

MA a B a B 

where the (nonvanishing) macrovariable /i' A is the average spacing of the energy levels, n' B and a' B are the average 
intensity and variance, respectively of the magnetic energy arising from the interaction of the tilted magnetic field 
with the spin | particle magnetic moment. In this case, I show that asymptotic behavior of £(^ haotlc ) j s lj ne ar in r 
(linear IGE growth), 

The equations for l 5(^ lte s rable ) anc j ^(chaotic) arg ^ e information-geometric analogue of the entanglement entropies 
defined in standard quantum information theory in the regular and chaotic cases, respectively. In addition, I emphasize 
that the statistical volume element Vm 3 (t) (see ©) may play a similar role as the computational complexity in 
conventional quantum information theory. These results warrant deeper analysis to be fully understood. 



IV. CONCLUSION 



In this paper I proposed a theoretical information-geometric framework suitable to characterize chaotic dynamical 
behavior of arbitrary complex systems on curved statistical manifolds. Specifically, an information-geometric analogue 
of the Zurek-Paz quantum chaos criterion of linear entropy growth and an information-geometric characterization of 
regular and chaotic quantum energy level statistics was presented. I hope that my work convincingly shows that this 
information-geometric approach may be useful in providing a unifying criterion of chaos of both classical and quantum 
varieties, thus deserving further research and developments. 

The descriptions of a classical chaotic system of arbitrary interacting degrees of freedom, deviations from Gaussianity 
and chaoticity arising from fluctuations of positively curved statistical manifolds are currently under investigation. I 
am also investigating the possibility to extend the IGAC to quantum Hilbert spaces constructed from classical curved 
statistical manifolds and I am considering the information-geometric macroscopic versions of the Henon-Heiles and 
Fermi-Pasta-Ulam /3-models to study chaotic geodesic flows on statistical manifolds. 
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